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We present branching fraction and CP asymmetry measurements as well as angular studies of 
B cf>(f>K decays using 464 x 10 6 BB events collected by the BABAR experiment. The branching 
fractions are measured in the <f><j> invariant mass range below the r/ c resonance (m^ < 2.85 GeV). We 
find B(B+ (f>(f>K + ) = (5.6±0.5±0.3) x 1(T 6 and B(B° (j><j>K°) = (4.5±0.8±0.3) x 10" 6 , where 
the first uncertaintiy is statistical and the second systematic. The measured direct CP asymmetries 
for the B ± decays are A C p = -0.10 ± 0.08 ± 0.02 below the n c threshold (m w < 2.85 GeV) and 
Acp = 0.09 ± 0.10 ± 0.02 in the n c resonance region (m</,0 in [2.94,3.02] GeV). Angular distributions 
are consistent with J p = 0~~ in the r\ c resonance region and favor J p — + below the n c resonance. 
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The violation of CP symmetry is a well-known require- 
ment for the matter-antimatter imbalance of the uni- 
verse The BABAR [1| and Belle Q experiments at the 
high-luminosity B factories, PEP-II 0] and KERB @, 
have made numerous CP asymmetry measurements us- 
ing datasets two orders of magnitude larger than their 
predecessors. All of these measurements are consistent 
with a single source of CP violation - the complex phase 
within the CKM quark mixing matrix of the Standard 
Model [6]. However, with the small amount of CP viola- 
tion from the CKM matrix, it is difficult to explain the 
matter-antimatter asymmetry of the universe [7|. This 
motivates searches for new sources of CP violation. 

A method to search for new sources of CP-violating 
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phases is to measure CP asymmetries in hadron decays 
that are forbidden at the tree level Q. Since the lead- 
ing decay amplitude is a one-loop process, contributions 
within the loop from virtual non-Standard-Model parti- 
cles cannot be excluded. The quark interactions with the 
non-Standard- Model particles can introduce new CP vi- 
olating phases in the decay amplitude, which can lead 
to observable non-zero CP asymmetries. Decays of B 
mesons with a b — ¥ sss transition have been extensively 
studied for this reason. 

The three body B — ► <fx\>K decay is a one-loop "pen- 
guin" b — > sss transition. This final state can also oc- 
cur through the tree-level decay B — > r\ c K, followed 
by r\ c — > 4>4>, where the B decay is a b — > ccs transi- 
tion. If the <jxj) invariant mass m§§ in the three-body 
B — > <f><f>K decay is close to the i] c resonance, the tree and 
penguin amplitudes may interfere. Within the Standard 
Model, the relative weak phase between these amplitudes 
is &rg(VtbV t * s /V c bV* s ) w 0, so no CP violation is expected 
from the interference. However, new physics contribu- 
tions to the penguin loop in the B — > cj)(j)K decay could 
introduce a non-zero relative CP violating phase, which 
may then produce a significant direct CP asymmetry Q . 
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Measurement of a significant, non-zero direct CP asym- 
metry would be an unambiguous sign of new physics. A 
previous measurement of the direct CP asymmetry [l(3] 
was consistent with zero, but was also limited by a large 
statistical uncertainty. The B + and B° branching frac- 
tions have been previously measured [l(| [H| to be a few 
times 10~ 6 . Theoretical predictions of the branching 
fractions are of the same order [12[ [13| ■ 



I. DATASET AND DETECTOR DESCRIPTION 

We present measurements of the B + — > cfxf>K + and 
B° — > (j)(j)K branching fractions [l4[ and direct CP 
asymmetry A C p = [N(B-)~N{B + )}/[N(B-) + N(B+)} 
as well as studies of angular distributions performed us- 
ing 464 x 10 6 BB pairs collected by the BABAR experiment 
at the SLAC National Accelerator Laboratory. The di- 
rect CP asymmetry is measured both below and within 
the rj c resonance region of the 4>4> invariant mass with 
these regions defined as < 2.85 GeV and within 
[2.94, 3.02] GeV, respectively p]|. The branching frac- 
tions are measured in the region below the rj c reso- 
nance (m00 < 2.85 GeV). 

The BABAR detector is described in detail elsewhere ■ 
What follows is a brief overview of the main features of 
the detector. The detector has a roughly cylindrical ge- 
ometry, with the axis along the beam direction. The tra- 
jectories, momenta, and production vertices of charged 
particles are reconstructed from position measurements 
made by a silicon vertex tracker (SVT) and a 40-layer 
drift chamber (DCH). The SVT consists of 5 layers of 
double-sided silicon strip detectors which provide preci- 
sion position measurements close to the beam interaction 
region. Both the SVT and DCH measure the specific 
energy loss (dE/dx) along the charged particle trajec- 
tory, which is used to infer the particle mass from the 
velocity dependence of the energy loss and the momen- 
tum measurement. The tracking system is inside a uni- 
form 1.5 T magnetic field provided by a superconduct- 
ing solenoid. Outside the tracking system, an array of 
quartz bars coupled with an array of phototubes (DIRC) 
detects the Cherenkov light produced when a charged 
particle travels through the quartz bars. The measured 
Cherenkov angle is used to infer the particle mass from 
the velocity dependence of the Cherenkov angle and the 
measured momentum. The energies of photons and elec- 
trons are determined from the measured light produced 
in electromagnetic showers inside a Csl crystal calorime- 
ter (EMC). Gaps in the iron of the magnet flux return 
are instrumented with resistive plate chambers and lim- 
ited streamer tubes, which are used to identify muons 
and neutral hadrons (IFR). 

We use Monte Carlo (MC) samples to determine the 
signal selection efficiency. The MC events are generated 
with EvtGen [lj| and simulated using Geant4 [l7j |. 



II. EVENT SELECTION 

We select events containing multiple hadrons by re- 
quiring at least three charged tracks in the event and the 
ratio of the second to zeroth Fox- Wolfram [l8j moments 
i?2 to be less than 0.98. 

Charged kaon candidates are required to pass a selec- 
tion based on a likelihood ratio which uses the SVT and 
DCH dE/dx and the DIRC Cherenkov angle measure- 
ments as inputs to the likelihood. The ratio is defined as 
Rhh' = £h/(£h + £h>), where h and b! are K, n, or p. 
The minimum kaon selection criterion is Rk-k > 0.2 or 
Rpn > 0.2. This selection has an efficiency greater than 
98% for kaons and a pion efficiency of less than 15% below 
a lab momentum of 2.5 GeV. Candidate <f> — > K + K~ de- 
cays are constructed from oppositely-charged kaon can- 
didates with an invariant mass in the range of 0.987 to 
1.2 GeV. At least one of the kaons in each — > K + K~ 
candidate must also satisfy the more stringent criteria of 
Rr-k > 0.8176 and Rk p > 0.018, which has an efficiency 
greater than 90% for kaons and a pion efficiency of less 
than 3% below a lab momentum of 2.5 GeV. 

Candidate Kg — > tt + it~ decays are constructed from 
oppositely-charged pion candidates with an invariant 
mass in the range of 0.486 to 0.510 GeV. The pion tracks 
are fit to a common vertex. The \ 2 probability of the 
vertex fit must be greater than 0.001. The typical ex- 
perimental resolution on the measured K§ flight length 
in the plane transverse to the beam is around 0.2 mm 
or less. We require the transverse flight length to be at 
least 2 mm. 

Candidate B — > cfxj)K decays are constructed from 
pairs of <f) candidates that do not share any daughters 
and either a K^ or a Kg candidate. The cj> and if ± 
candidates are constrained to a common vertex. We re- 
ject combinatoric background by requiring the B can- 
didate to have kinematics consistent with T(AS) —> BB 
using two standard variables: toes and AE. The energy- 
substituted B mass is defined as toes = V^beam — "P* 2 1 
where E£ cam andp* are the beam energy (yfs/2) and the 
reconstructed B momentum, both in the center-of-mass 
(CM) reference frame. The energy difference is defined 
as AE = E* — E^ cam , where E* is the reconstructed B 
energy in the CM frame. We require toes and AE to be 
within [5.20, 5.29] GeV and [-0.1,0.1] GeV, respectively. 
The experimental resolution is about 2.7 MeV for Toes 
and 15 MeV for AE. 

The toes interval includes a large "sideband" region 
below the area where the signal is concentrated near the 
B mass. The AE interval also is wide enough to include 
sideband regions where the signal probability is very low. 
Including events in the sideband regions enables us to 
determine the probability density functions (PDFs) of 
the combinatoric background directly in the maximum 
likelihood (ML) fits of the data. 

About 7% of events in signal Monte Carlo samples have 
more than one B — > <f>(f>K candidate. If there are multi- 
ple B + —> (f>(f>K + candidates in a single event, we select 



7 



the B + candidate with the smallest mass Xm defined as 

53 i ( ~ m '~"'°" ) , where the sum i is over the two 4> can- 
didates, m,i (too) is the reconstructed (nominal) (f> mass, 
and a is the RMS of the reconstructed distribution for 
properly reconstructed <f> candidates. If there are more 
than one B + candidates that use the same two <fi can- 
didates, we choose the B + candidate with the highest 
quality K + identification for the K + from the B + de- 
cay. If the quality level of the K + identification is the 
same for these B + candidates, we choose the B + candi- 
date with the highest vertex \ 2 probability. For events 
with multiple B° — > c/xpKg, the sum for Xm includes the 
Kg — > 7r + 7r~, and we choose the B° with the smallest 
Xm- F° r both the B + and B° decay modes, the proba- 
bility that the the algorithms described above choose the 
correct candidate is about 87%. 

The reconstruction and selection efficiencies for events 
with m^tf, < 2.85 GeV are determined from Monte Carlo 
samples to be 28.0% and 22.5% for the B+ and B° modes, 
respectively. 

We use control samples of B — > D S D decays where 
D s — > 4>tt, D — > Ktt, and <fi — > K + K~ to determine 
corrections to the B — > (jxpK signal probability density 
function parameters determined from Monte Carlo sam- 
ples in the maximum likelihood fits described below. 

A. Continuum Background 

The events that pass the selection above with at 
least one B candidate are primarily background events 
from the continuum (e + e _ — > qq with q = u,d,s,c). 
We reduce this background by using a Fisher discrim- 
inant (J 7 ), which is the linear combination of seven 
variables and is optimized for maximum separation 
power of signal and the continuum background. The 
seven variables are listed below. These variables are 
commonly used by the BABAR experiment in analyses 
of charmless B decays, where the primary background 
is from continuum events. They take advantage 
of aspects of the production distributions and event 
topologies of BB versus continuum qq production events. 

• | At / crat | : the absolute value of the reconstructed 
proper time difference between the two B decays 
divided by its uncertainty [19(. 

• |-FT|: the absolute value of the standard BABAR 
flavor tagging neural network output [l9| . 

• | cos 0j* h | : the absolute value of the cosine of the 
angle between the B candidate thrust axis and the 
thrust axis of the rest of the event computed in the 
CM frame. The thrust axis is the direction that 
maximizes the scalar sum of the projection of the 
track momenta on that direction. 

• | cos 9 Bthr | : the absolute value of the cosine be- 
tween the thrust axis of the B candidate and the 




Fisher discriminant 

FIG. 1: Distributions of the Fisher discriminant T for B — > 
(fxf> K signal (solid black) and e + e~ — > qq continuum (dashed 
red) Monte Carlo simulation. 



beam axis in the CM frame. Signal events have a 
uniform distribution in this variable, while contin- 
uum background follows a |1 + cos 2 0\ distribution, 
where 9 is the angle between the thrust direction 
and the beam axis. 

• | cos B | : the absolute value of the cosine of the 
angle between the B direction and the beam axis 
in the CM frame. The angular distribution of the 
signal follows a sin 2 8* B distribution, while the con- 
tinuum background is uniformly distributed. 

• Lo and L?: The zeroth and second angular mo- 
ments of the momentum flow of the rest of the 
event about the B thrust axis, defined as Lj = 
^2nPi\ cosfjp, where the angle 9i is the angle be- 
tween track i and the B thrust axis and the sum 
excludes the daughters of the B candidate. The 
calculations are done in the CM frame. 

Distributions of T for signal and continuum MC samples 
are shown in Fig. [TJ The Fisher discriminant T is used 
as one of several variables in the maximum likelihood fits 
described below. 



B. Peaking Backgrounds 

The ultimate detected state of our B — > <jxj>K sig- 
nal decay is five kaons. In addition to the <p reso- 
nance, there may be contributions to each K + K~ pair 
either from other intermediate K + K~ resonances, such 
as the Jo (980), or from non-resonant K + K~ contribu- 
tions. We use the K + K~ mass sidebands for each (j) 
candidate to determine the amount of B mesons that 
decay to the detected five-kaon state (which we denote 
B —> 5K) that are not coming from B —> <f>cj)K. The 
specific B decays that we consider as backgrounds are 
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B -> cj)K+K-K, B -> K+K-K+K-K, B -> f <j>K, 
and B — > f K + K~ K . The branching fractions for these 
decays are currently unknown. We call these B decays 
"peaking backgrounds" because properly reconstructed 
B candidates are indistinguishable from our B — > (fxpK 
signal in the toes, AE, and T variables. 

We perform unbinned extended maximum likelihood 
fits to determine the signal and combinatoric background 
yields and, in some cases, the charge asymmetry. All 
of the fits use the product of one-dimensional PDFs of 
toes, AE, and T in the likelihood. For the B — > 4><j>K 
branching fraction measurements, we also include PDFs 
for the invariant mass of each (j> — > K + K~ candidate 
(m^i and m^ 2 )- 

As a first step, we divide the m^i vs. to^2 plane [20j 
in the range of 0.987 to 1.200 GeV into five mutu- 
ally exclusive zones. We fit for the B — > 5K yield in 
each zone using only toes, AE, and T in the likeli- 
hood. The zones are based on various combinations of 
the 4> signal and sideband regions, which are defined as: 
Low-SB [0.987,1.000] GeV, phi-signal [1.00,1.04] GeV, 
and High-SB [1.04,1.20] GeV. Each of the five zones is 
chosen so that either the B — > (jxpK signal or one of the 
four peaking B backgrounds is concentrated in the re- 
gion. We compute the number of peaking background 
events within the to^ range used for the branching frac- 
tion fit by using the results of the five zone fits as de- 
scribed below. 

Figure [2] shows the distribution of events in the 
vs. 77i0i plane for the selected B + — > 5K candidates in 
the data. To enhance the B + — > 5K signal for the fig- 
ure, we have required 7n E s > 5.27 GeV, \AE\ < 0.040 
GeV, and T < 0.0. The inset of the figure shows the 
definition of the five zones. A concentration of events in 
the phi-signal region for both <f) candidates (zone 1) is 
clearly evident. The region defined as phi-signal com- 
bined with High-SB for either <f> candidate (zone 2) con- 
tains the largest fraction of the B — > <{)K + K~K mode, 
although the B — > (jxpK signal also populates this region 
due to cases where one <fi is mis-reconstructed. The zone 
where the invariant mass of both <j) candidates is in the 
High-SB region (zone 3) contains the largest concentra- 
tion of the non-resonant B — > K + K~ K + K~ K mode. 
Zones 4 and 5 contain a large fraction of the B — > fo4>K 
and B —> f$K + K~ K modes, respectively, and very small 
fractions of the other three modes. 

Monte Carlo samples for the five B decay modes (sig- 
nal plus four peaking background modes) are used to 
determine the fraction of events in each zone (i) for each 
decay mode (J), which we denote with the matrix fy. 
The total B 5K yield (rij) is determined for each zone 
i using five separate maximum likelihood fits of the data. 
The yield for each B decay mode (Nj) and the amount 
of each mode j in zone i (ijy) can be determined from 

Nj = fij 1 ™* and n v = fii N 3- C 1 ) 

i 

Zone 1 corresponds to the to^ range used in the branching 
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FIG. 2: Data distribution of B + — > 5K events in the vs. 
m$\ plane for events with m^^ < 2.85 GeV. To enhance the 
B — > 5K signal for the figure, we have required rriES > 5.27 
GeV, \AE\ < 0.040 GeV, and T < 0.0. The efficiency of these 
additional requirements, relative to the nominal selection, is 
about 70% for the signal. The inset shows the definition of 
the five zones. 



fraction maximum likelihood fit. 



C. Maximum Likelihood Fits 

The extended maximum likelihood fits in the five zones 
determine the B — ► 5K signal and combinatoric back- 
ground yields in each zone. The B — > 5K signal is split 
into properly reconstructed and misreconstructed ( "self- 
crossfeed") components, with the self-crossfeed fraction 
fixed. The self-crossfeed component is defined as events 
where a true B — > 5K decay is present in the event, but 
one or more tracks used in the reconstructed B are either 
from the other B in the event or not real. In zone 1, the 
self-crossfeed fraction for B — > <jxf>K decays is around 7%. 

The properly reconstructed B — > 5K signal component 
is described by the following PDFs: a Crystal Ball func- 
tion [2l| for toes, the sum of three Gaussians for AE, 
and the sum of a bifurcated Gaussian and a Gaussian 
for T . The Crystal Ball function is a Gaussian modi- 
fied to have an extended power-law tail on the low side. 
The B — > 5K signal PDF parameters are determined 
from MC samples with corrections to the Toes and AE 
core mean and width parameters from the B — > D S D 
control samples. The mean corrections are 0.04 ±0.11 
MeV and —3.5 ±0.8 MeV for toes and AE, respectively. 
The width scale factors are 1.10 ± 0.04 and 1.04 ± 0.05 
for Toes and AE, respectively. The combinatoric back- 
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ground is described by the following PDFs: an empirical 
threshold function [22j for ttiesj a first-order polynomial 
for AE, and the sum of two Gaussians for T. Most of 
the combinatoric background PDF shape parameters are 
determined in the fits. 

The results of the five zone fits for the B + and B° 
modes are given in Tables [V] and IVI| respectively, in the 
appendix. The B — > (fxfiK signal is observed in both 
the B + and B° samples. The B° —> (fxfiK decay has 
not been observed previously. The B yield in zone 2 for 
the B + mode is significant, but about half of this is due 
to misreconstructed B — » <ficj)K signal. The computed 
B+ (f)K+K-K+ and B+ -> K+K-K+K-R+ yields 
are positive, but the significance is less than two standard 
deviations. There is no evidence of either B — > fo4>K 
or B — > foK + K~ K . The branching fraction maximum 
likelihood fits use the range that corresponds to zone 
1. We fix the yield of each of the four peaking back- 
ground modes to the zone 1 value in Table [V] or I VII for 
the branching fraction fit described below. 



III. BRANCHING FRACTION ANALYSIS 

The maximum likelihood fit used to measure the B — > 
(j)cj)K yield below the rj c resonance for the branching 
fraction measurement restricts the event selection with 
m<p<p < 2.85 GeV and within [1.00,1.04] GeV, which 
corresponds to zone 1 in the peaking background discus- 
sion above. The fit components are B —> cfxpK signal, 
combinatoric background, and the four peaking back- 
grounds. 

In addition to toes: AE, and T , PDFs for m^i and 
are included in the likelihood function. For each fit 
component, each <fi candidate has a PDF that is the sum 
of a properly reconstructed 4> — > K + K~ decay, given 
by a relativistic Breit-Wigner function, and a misrecon- 
structed </>, described by a first-order polynomial. The 
m^i and ra^i PDFs are combined in a way that is sym- 
metric under 1 2 exchange and takes into account the 
fractions of events where both </> candidates are prop- 
erly reconstructed, one </> is misreconstructed, and both 
<p candidates are misreconstructed. 

In addition to the signal and combinatoric background 
yields, the charge asymmetry for the signal and combi- 
natoric background components and most of the combi- 
natoric background PDF parameters are determined in 
the fit. 

The results of the B + and B° fits are shown in Figs. [3] 
and 21 respectively. To reduce the combinatoric back- 
ground in each distribution shown in the figures, a re- 
quirement is made on a likelihood ratio, which is based 
on all the fit variables except the one plotted. 

The fitted charge asymmetry Acp for the background 
component is 0.02 ± 0.03. The charge asymmetry for the 
signal component is —0.10 ± 0.08. The fitted yields of 
B + and B° signal candidates with < 2.85 GeV are 
178 ± 15 events and 40 ±7 events, respectively, where the 




1 2 
Fisher 



FIG. 3: Results of fitting the B + -> 4>(t>K + sample for 
< 2.85 GeV. The dashed red curve is the sum of the 
combinatoric and peaking background components. The solid 
blue curve is for all components. A requirement on a likeli- 
hood ratio based on all fit variables except the one plotted is 
made to reject most of the background. The likelihood ratio 
requirements are about 84% efficient for the signal. 




FIG. 4: Results of fitting the 73° (j><f)K% sample for 
< 2.85 GeV. The dashed red curve is the sum of the 
combinatoric and peaking background components. The solid 
blue curve is for all components. A requirement on a likeli- 
hood ratio based on all fit variables except the one plotted is 
made to reject most of the background. The likelihood ratio 
requirements are about 84% efficient for the signal. 



uncertainties are statistical only. 
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A. Systematic Uncertainties 



Table U summarizes the systematic uncertainties on the 
B —> <fi(f>K branching fractions in the < 2.85 GeV 
region. The systematics are divided into additive uncer- 
tainties that affect the B yield measurement and multi- 
plicative uncertainties in the branching fraction calcula- 
tion. 

The uncertainties from the corrections applied to the 
PDF parameters such as the tties and AE core mean and 
width for the signal component, which are derived from 
data control samples, are listed under "ML Fit Yield". 
The signal Fisher and m<f, core Gaussian mean and width 
parameters are not corrected in the fit, because data con- 
trol sample measurements are consistent with the Monte 
Carlo. However, we did vary the signal Fisher and 
core Gaussian mean and width parameters by the sta- 
tistical uncertainty of the data control sample measure- 
ments. These variations are also included under "ML 
Fit Yield" . The fit bias systematic is taken to be half of 
the bias correction added in quadrature with the statis- 
tical uncertainty on the bias. We vary the fixed peaking 
background yields by their statistical uncertainties (see 
Tables [V] and IVij) and by varying the fractions /„■ . The 
fixed self-crossfeed fraction for the signal component was 
varied by ±2%. Adding the individual uncertainties in 
quadrature, the total additive systematic uncertainties 
on the B + and B° signal yields are 6.2 and 1.8 events, 
respectively. 

The uncertainty on the track reconstruction efficiency 
is ±0.23% per track, which is taken to be fully corre- 
lated for the charged kaons. The Kg reconstruction ef- 
ficiency has an uncertainty of 1.5%. The 4> — » K + K~ 
and K% — > 7r + 7r~ branching fractions are taken from the 
PDG [23| and are varied by their one standard devia- 
tion uncertainties. The systematic uncertainty on the 
K^ identification criteria was estimated by comparing 
the ratio of the B yield with the nominal selection to 
the B yield requiring all K^ to pass the tighter selection 
in the data and the MC samples. This gives an uncer- 
tainty of 3% for the B + mode and 2% for the B° mode. 
Adding the individual uncertainties in quadrature, the 
overall multiplicative systematic uncertainties are 3.6% 
for the B + mode and 3.2% for the B° mode. 

The signal charge asymmetry has been corrected for 
a bias due to differences in the K + and K~ efficiencies 
by adding +0.010 ± 0.005 to the asymmetry. The overall 
2% systematic uncertainty takes into account uncertain- 
ties on the charge dependence of the tracking efficiency, 
material interaction cross section for kaons, and particle 
identification. 



TABLE I: Summary of the systematic and statistical uncer- 
tainties for the branching fraction measurements. 



Quantity B + - 


-> cj>cj>K + B { ' 


-> 4>4>K" 


Fit Stat. Uncertainty (events) 


15.1 


7.0 


Additive Uncertainties 


(events) 




ML Fit Yield 


3.3 


1.0 


ML Fit Bias 


1.6 


0.2 


Peaking BG, region fits 


3.5 


1.2 


Peaking BG, fij values 


3.0 


0.8 


Self-Crossfeed Fraction 


1.8 


0.4 


Total Additive Syst. (events) 


6.2 


1.8 


Multiplicative Uncertainties (%) 




Tracking Efficiency 


1.2 


1.0 


K® Reconstruction Efficiency 




1.5 


Number BB 


1.1 


1.1 


B K+K~) 


1.2 


1.2 


B {K° s -> tt+tT) 




0.1 


MC Statistics 


0.1 


0.1 


m^rj, Cut Efficiency 


0.2 


0.3 


K Identification 


3.0 


2.0 


Total Multiplicative Syst. (%) 


3.6 


3.2 


Total Systematic [B] (xlCT b ) 


0.3 


0.3 


Statistical [B] (xlO -6 ) 


0.5 


0.8 



B. Branching Fraction Results 

Table [TT] summarizes the B — > (fxfiK branching fraction 
results for < 2.85 GeV. We find 

B(B+ -> cjxj)K + ) = (5.6 ±0.5 ±0.3) x 10~ 6 
B{B° -> 4>4>K°) = (4.5 ±0.8 ±0.3) x 10~ 6 , 

where the first uncertainty is statistical and the second 
systematic. These results are consistent with and super- 
sede the previous measurements [ll[ by the BaBAR Col- 
laboration. The Belle collaboration measurements [13] 
are lower, though they are statistically compatible. Our 
branching fraction measurements are higher than the the- 
oretical predictions of [l2| and [l3T |. 

IV. CP ASYMMETRY IN rj c RESONANCE 
REGION 

As was mentioned in the introduction, a significant 
non-zero direct CP asymmetry in the r\ c resonance region 
of m^^ would be a clear sign of physics beyond the Stan- 
dard Model. For this measurement, we use the simpler 
likelihood, based on m^s, AE, and T. Figure [5] shows 
the fitted B + — > (j><f>K + yield as a function of m^. The 
r\ c resonance is clearly visible. Narrow bins around the 
Xco and \c2 resonances do not show a significant excess 
above the broad non-resonant component. 

The results of fitting the events in the range of 
[2.94,3.02] GeV are given in Table HID For the CP asym- 
metry, we find 

Acpijn^ in [2.94,3.02] GeV) = -0.09 ± 0.10 ± 0.02, 
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TABLE II: Branching fraction and charge asymmetry results 
for B — > (jxj>K in the region < 2.85 GeV. The statisti- 

cal significance is given by ln(£ max /£o), where £ max is 
the maximum likelihood and Co is the likelihood for the hy- 
pothesis of no (jxj>K signal. The significance does not include 
systematic uncertainties. 





B + -> 4>4>K + 


B" -> cjxj>K u 


Events to fit 


1535 


293 


Fit signal yield 


178 ± 15 


40 ± 7 


ML-fit bias (events) 


3.0 ± 0.5 


0.0 ± 0.2 


MC efficiency (%) 


28.0 


22.5 


IWi(%) 


24.2 


8.4 


Stat, significance 


24 


11 


/3(kt 6 ) 


5.6 ± 0.5 ± 0.3 


4.5 ± 0.8 ± 0.3 



Signal Acp 
Comb. Bkg. Acp 



-0.10 ±0.08 ±0.02 
0.02 ± 0.03 
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FIG. 5: Fitted _B + — » <fxf>K + yield as a function of m,^. 
Each point shows the results of a maximum likelihood fit of 
the events in that bin. The inset is the same data with an 
expanded vertical range to show the shape of the non-resonant 
component more clearly. The yield has been divided by the 
bin width and scaled by 0.027 GeV, which is the bin width 
of the three bins in the r\ c resonance region ([2.94,3.02] GeV 
and dashed vertical lines in the inset). The two narrow bins 
above the r) c are centered on the XcO (bin range [3.400,3.430] 
GeV) and the \c2 (bin range [3.552,3.560] GeV). 



TABLE III: Fit results for B + — > <jxj>K + within r\ c resonance 
region (m^ within [2.94,3.02] GeV). The signal charge asym- 
metry Acp has been corrected by adding ±0.010 ± 0.005 to 
the fitted asymmetry. 



ML fit quantity /Analysis B + -» (fxj)K + 


Events to fit 


181 


Fit signal yield 


100 ± 10 


MC efficiency (%) 


29.2 


Corr. Signal Acp 


-0.09 ±0.10 ±0.02 


Comb. Bkg. Acp 


-0.06 ±0.11 



would expect 62 ± 22 signal events, ignoring the non- 
resonant B + — > <{xj)K + contribution and any interfer- 
ence between the resonant r\ c and non-resonant ampli- 
tudes. We do not use our B + event yield to measure 
B(B + i] c K + ) x B(r] c -t (fx/)) due to the potentially 
large interference effects between the resonant and non- 
resonant (jxj) amplitudes which we can not easily quantify. 

The Acp may integrate to zero, even if there is a con- 
tributing non-Standard-Model amplitude with a non-zero 
CP violating phase. However, in this case the phase vari- 
ation of the r) c resonance amplitude could give non-zero 
Ac p values with opposite signs above and below the peak 
of the resonance. We have performed the measurement in 
two ranges, splitting the r/ c region into two regions (above 
and below the peak of the resonance) . The results are 



A C p{m u in [2.94,2.98] GeV) 
Acpim^ in [2.98,3.02] GeV) 



-0.10 ±0.15 ±0.02 
-0.08 ±0.14 ±0.02, 



both of which are consistent with zero, as expected in the 
Standard Model. 



V. ANGULAR STUDIES 

We use the angular variables that describe the B + — > 
4>cf>K + decay to investigate the spin components of the 
4>(f> system below and within the rj c resonance. The angles 
are defined as follows. 



where the first uncertainty is statistical and the second 
uncertainty is systematic. The value above includes the 
same 1% bias correction and has the same 2% overall 
systematic uncertainty as the signal charge asymmetry 
below the r\ c resonance as described above. 

The fit yields 100 ± 10 signal candidates. Using 
B(B+ -> T) C K+) = (9.1 ± 1.3) x 10" 4 and B(rj c -> 00) = 
(2.7 ± 0.9) x 10~ 3 from the PDG @, a B + 
4> — > K + K~ reconstruction efficiency of 29% in the 
r\ c resonance region, and an efficiency of 78% for the 
77i00 window of [2.94,3.02] GeV for the rj c resonance, we 



• Gi, (i = 1, 2) : The 0i angle is the angle between the 
momentum of the K + coming from the decay of <j>i 
in the <pi rest frame with respect to the boost di- 
rection from the (jxj) rest frame to the (f>i rest frame. 

• x : The x angle is the dihedral angle between the 
4>i and (f>2 decay planes in the (jxj) rest frame. 

• 000: The angle is the angle between one of 
the (j) mesons in (jxj) rest frame with respect to the 
boost direction from the B + rest frame to the (jxj) 
rest frame. 
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TABLE IV: Quality of the angular fits shown in Fig. [7] The 
first column is the interval for the events in the fit. The 
last column is the p-value of the \ 2 goodness-of-fit test for 
the hypothesis indicated in the third column. 



a 100 



(GeV) 


Variable 


rUr 


1 AT 

X /Ndol 


X prob. 


To O/i o fiOl 

[z.y4, o.Uzj 


X 


sin'' x 


y.ol/y 


U.oy 


[2.94, 3.02] 


X 


uniform 


60.3/9 


1.2 x 10~ 9 


< 2.85 


X 


sin 2 x 


41.6/9 


3.9 x 10 -6 


< 2.85 


X 


uniform 


18.5/9 


0.030 


[2.94, 3.02] 


cos 9 


sin^x 


9.97/9 


0.39 


[2.94, 3.02] 


cos 9 


uniform 


60.5/9 


1.1 x 10~ 9 


< 2.85 


cos 9 


sin 2 x 


32.9/9 


1.7 x 10~ 4 


< 2.85 


cos 9 


uniform 


25.8/9 


2.2 x 10" 3 


[2.94, 3.02] 


cos Ofal 


uniform 


9.02/9 


0.44 


< 2.85 


cos 


uniform 


5.01/9 


0.83 



We project the J p = component by making a his 
togram of TO00 weighting each event by 

P 2 (cos6 1 ) Re[Y 2 2 (9 2 , X )] = 
— {3 cos 6\ — 1) sin 9 2 cos 2%, 



(2) 



where P 2 is a second-degree Legendre polynomial and 
Y 2 2 is a spherical harmonic with £ = 2 and m — 2. In 
each bin, the J p = 0~ component yield is projected 
out, while the combinatoric background averages to zero. 
To do this, we select events in a signal region defined 
by: m ES > 5.27 GeV, |A.E| < 40 MeV, m within 
[1.01,1.03] GeV, and T < 0.5. The efficiency of these re- 
quirements, relative to the selection used in the asymme- 
try measurement, is about 78% for signal events and 2.9% 
for combinatoric background. The combinatoric back- 
ground that remains after this selection is shown using 
data events in the sideband region (toes < 5.27 GeV and 
\AE\ < 100 MeV) scaled by 0.065, which is the signal- 
to-sideband ratio for the combinatoric background. 

The results are shown in Fig. [TJ] The weighted yield in 
the 7y c region is consistent with all of the B + — > <jxj)K + 
events having J p — Q~. Just below the rj c region, the 
weighted yield is consistent with zero. The excess in the 
bins near 2.2 GeV may be due to the 7?(2225) seen in 
J/ip -t 700 events at Mark III [13] and BES 

Figure [7] shows background-subtracted distributions of 
X, cos#i, and [cos 6*00 | for the nominal event selection. 
The background subtraction is done with the technique 
described in reference (2(|. Since there is no meaningful 
distinction between 0i and <p 2 , we combine the cos 9\ and 
cos 9 2 distributions into one plot of cos 9. The reconstruc- 
tion and selection efficiency, determined from MC sam- 
ples, is uniform in x and cos#i, but not in | cos #00 1, so 
the | cos 6*00 1 distribution is efficiency corrected. For each 
distribution, we performed a simple least-x 2 fit to the dis- 
tributions expected for both J p = 0~ and J p = + for 
the 00 system. 

For a J p — 0~ state, we expect x to have a sin 2 x 
distribution, while x should be uniform for J p = + . 
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FIG. 6: Histograms (top) and weighted distributions (bot- 
tom) of for the signal region (solid points) and data 
sideband selection (dashed with open diamonds) defined in 
the text. The sideband distributions have been normalized 
to the expected level of combinatoric background remaining 
after the signal region selection. Events in the bottom dis- 
tribution were weig hted by P 2 (cos6>i) Re [Y 2 2 (0 2 ,x)] which 
projects out the J p = 0~ component. The yield has been 
divided by the bin width and scaled by 0.21f GeV, which is 
the width of the bin covering the r/ c resonance ([2.875,3.086] 
GeV). 



The signal events in the rj c resonance region are consis- 
tent with a sin 2 x distribution while the signal below the 
i] c resonance is not. For a J p = 0~ state, the distribu- 
tions of cos 6i are expected to have sin 2 9i distributions, 
while a J p = + state is expected to have uniform cos 9{ 
distributions. The events in the r\ c resonance region are 
consistent with a sin 2 9i distribution, while the events 
below the r\ c resonance show a deviation from a sin 2 9i 
shape. 

Finally, a spin-zero state should have a uniform 
| cos 6*00 1 distribution. The efficiency-corrected distribu- 
tions shown in Fig. [7j both within and below the T) c res- 
onance region, are consistent with a uniform | cos 9^\ 
distribution. 
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m # in [2.94,3.02] GeV m # in [2.94,3.02] GeV nu, in [2.94,3.02] GeV 




FIG. 7: Background subtracted angular distributions in the rj c resonance region (m^ in [2.94,3.02] GeV for the top row) and 
below the rj c resonance (m^ < 2.85 GeV for the bottom row). The reconstruction and selection efficiency is uniform in \ 
(left) and cos 8 (center), but dependent on | cos#^| (right), so the right column has been efficiency corrected. The red dashed 
line shows a least-^ 2 fit of the points to a uniform distribution while the solid black curve shows a fit to the expectation for a 
J p — 0~ state decaying to <jxj). 



VI. SUMMARY AND CONCLUSIONS 

We have measured the branching fractions and charge 
asymmetries of B — > cjxfiK decays below the rj c resonance 
in the (jxj) invariant mass (m^ < 2.85 GeV). We observe 
both B + — > (jxj)K + and B° -> <t><l>Kg, each with a sig- 
nificance of greater than five standard deviations. The 
B° —> 4><fiKg decay has not been observed previously. 
Our branching fraction measurements are higher than 
the theoretical predictions of [l2[ and fl3j . 

We have measured the charge asymmetry for B + — > 
(j)(f)K + in the r/ c resonance region, where a significant 
non-zero value would be an unambiguous indication of 
new physics. Our measurement is consistent with zero, 
which is the expectation of the Standard Model. 

Finally, we have studied the angular distributions of 
B + — > (fxf>K + decays below and within the r\ c resonance. 
We conclude from these studies that the non-resonant 
B + — > <f>4>K + events below the r\ c resonance are, on aver- 
age, more consistent with J p = + than J p = CP, while 
the distributions within the r\ c resonance region are all 
consistent with J p = 0~. 
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Appendix A: Peaking background zone fit results 

The results of the fits of the B — > 5K yield in the five 
zones in the vs m^i plane (shown in Fig. [2]) and 
the derived yields for each of the five B decay modes are 
given in Tables IVl and I VII below. 
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TABLE V: B+ -> 5K yield fit results for the five zones (B + yield column) and yields derived from the fraction matrix fij. 
The last row is computed from the inverted fraction matrix f~^ and the B + yield column. The remaining 5x5 yield matrix 
was computed from the last row and the fraction matrix fij . The uncertainties are from propagating the statistical errors from 
the fitted region yields (B + yield column) without including any uncertainties on the fraction matrix. 



Zone 


B + yield 




4>KKK ± 


KKKKK ± 




foKKK* 


1 


188.4 ± 16.0 


183.2 ± 17.0 


6.3 ±3.9 


0.8 ±0.8 


0±3.0 


0.3 ± 1.5 


2 


84.4 ± 18.0 


36 


39 


9 


-1.4 


2 


3 


49.7 ± 19.0 


3.8 


18 


26 


-0.2 


2.2 


4 


1.0 ±2.0 


1.3 


0.9 


0.2 


-1.7 


0.3 


5 


3.5 ±5.0 


0.2 


0.8 


1.3 


-0.2 


1.3 


1 - 5 




225 ±21 


65 ±40 


38 ±38 


-5.7 ±7.6 


6.2 ±26 



TABLE VI: B° — > 5K yield fit results for the five zones (B° yield column) and yields derived from the fraction matrix fij. 
The last row is computed from the inverted fraction matrix f~. x and the B° yield column. The remaining 5x5 yield matrix 
was computed from the last row and the fraction matrix fij . The uncertainties are from propagating the statistical errors from 
the fitted region yields (B° yield column) without including any uncertainties on the fraction matrix. 



Zone 


B u yield 




4>KKK" 


KKKKK" 




f KKK u s 


1 


43.4 ±8.0 


42.8 ±8.4 


1.2 ± 1.5 


0.1 ±0.3 


0±0.7 


0±0.02 


2 


15.1 ±7.0 


8 


7 


0.9 


-0.4 


-0.4 


3 


6.0 ±6.0 


1 


3 


2.5 


-0.1 


-0.4 


4 


0±0.5 


0.3 


0.2 





-0.5 





5 


0±0.8 





0.1 


0.1 





-0.2 


1 - 5 




52 ± 10 


12 ±15 


3.6 ± 12 


-1.6 ± 1.8 


-1.1 ±4.3 
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